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Abstract 

Centrifugal compressors are compatible with the low exit corrected flows found in the high pressure 
compressor of turboshaft engines and may play an increasing role in turbofan engines as engine overall 
pressure ratios increase. Centrifugal compressor stages are difficult to model accurately with RANS CFD 
solvers. A computational study of the CC3 centrifugal impeller in its vaneless diffuser configuration was 
undertaken as part of an effort to understand potential causes of RANS CFD mis-prediction in these types 
of geometries. Three steady, periodic cases of the impeller and diffuser were modeled using the TURBO 
Parallel Version 4 code: (1) a k-e turbulence model computation on a 6.8 million point grid using wall 
functions, (2) a k-e turbulence model computation on a 14 million point grid integrating to the wall, and 
(3) a k-co turbulence model computation on the 14 million point grid integrating to the wall. It was found 
that all three cases compared favorably to data from inlet to impeller trailing edge, but the k-e and k-co 
computations had disparate results beyond the trailing edge and into the vaneless diffuser. A large region 
of reversed flow was observed in the k-e computations which extended from 70 to 100 percent span at the 
exit rating plane, whereas the k-co computation had reversed flow from 95 to 1 00 percent span. Compared 
to experimental data at near-peak-efficiency, the reversed flow region in the k-e case resulted in an under- 
prediction in adiabatic efficiency of 8.3 points, whereas the k-co case was 1.2 points lower in efficiency. 

Nomenclature 


CFD 

computational fluid dynamics 

CMOTT 

Center for Modeling of Turbulence and Transition 

k 

turbulence kinetic energy (TKE) 

R 

radius 

Rte 

radius of impeller trailing edge (215.5 mm) 

RANS 

Reynolds-averaged Navier Stokes equations 

+ 

y 

dimensionless wall distance 

8 

dissipation of turbulence kinetic energy 

CO 

specific turbulence dissipation 
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Introduction 


Centrifugal compressors are frequently used in the high pressure compressor of turboshaft engines 
due to the low exit corrected flow rates in these machines. With the decreasing size of the high pressure 
compressor in turbofan engines as engine overall pressure increases, exit corrected flow rates may be 
sufficiently low such that the use of a centrifugal stage at the compressor exit may provide performance 
benefits over a typical axial compressor stage (Ref. 1). Thus, the ability to accurately simulate the flow 
within centrifugal stages is becoming increasingly important for compressor designers and researchers. 
However, Reynolds-averaged Navier-Stokes (RANS)-based computational fluid dynamics (CFD) results 
published in the literature tend to mis-predict the performance of centrifugal compressors when compared 
to available data (Refs. 2 to 5). In general, the computations tend to over-predict total pressure ratio, 
efficiency, and in the case of vaned diffuser configurations, the choking flow rate. 

The current research effort was undertaken in an attempt to investigate and understand the potential 
causes of the mis-prediction of centrifugal compressor impeller performance by RANS solvers. Figure 1 
shows the impeller studied in this investigation, a centrifugal compressor known as CC3. This test 
compressor is a 4.54 kg/s (10 lb m /s) machine which was scaled up from a 1.658 kg/s (3.655 lb m /s) Allison 
Engine design. The design exit corrected flow rate of the test compressor is 1.36 kg/s (3.0 lb m /s). The 
complete design description and hot shape geometry of the test compressor are given by McKain and 
Holbrook (Ref. 6). The test compressor was used in aerodynamics studies at NASA Glenn Research 
Center (Refs. 7 to 9). The impeller was tested in both vaned diffuser and vaneless diffuser configurations, 
and detailed measurements of the flow field were obtained. The data available for this test compressor 
qualifies it as a good candidate for further computational studies. 

The CFD solver used in the current work is the TURBO Parallel Version 4 code. TURBO is a three- 
dimensional, viscous, unsteady RANS solver for turbomachinery flows. Flux vector splitting is used in 
evaluating flux Jacobians on the left hand side, while Roe’s flux difference splitting is used to form a 
higher order TVD scheme for convective fluxes on the right hand side. The code can use multiple 
structured blocks with matching boundaries. Local time stepping can be employed to accelerate 
convergence for steady-state flows. Newton sub-iterations are applied to unsteady flows to obtain a 
converged solution within each time step (Ref. 10). Multiple blade rows can be modeled with sector 
periodic blade row meshes using a simple periodic boundary condition or with a single blade passage per 
blade row using phase lag boundary conditions (Refs. 1 1 and 12). The default turbulence model in 
TURBO is the NASA/CMOTT k-s turbulence model with wall functions (Ref. 13). The code has 
previously been validated on various geometries (Refs. 14 and 15). 



Figure 1. — CC3 impeller and vaned diffuser hardware (Ref. 9). 
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Figure 2. — Computed and experimental total-to-static 
pressure ratios of CC3 vaned diffuser configuration 
at design speed. 


As a precursor to the current research effort, unsteady phase-lag simulations of the CC3 impeller with 
the vaned diffuser configuration were generated using TURBO. The resulting speed line at the design 
speed is shown in Figure 2 alongside a comparison to the data. An over-prediction in choking flow rate 
was observed in the CFD result, and this speed line serves to illustrate a trend seen throughout the 
literature (Refs. 2 to 5). 

The CFD results shown are converged, unsteady phase-lag computations which were time-averaged 
over one revolution. Computed static pressures at the exit rating station were area-averaged. In these 
simulations, a grid of 4.3 million points with resolution deemed to be typical of turbomachinery CFD was 
used. A row of three support struts upstream of the impeller, as well as leakage flows, were neglected in 
the computation. The computation over-predicted pressure ratio by 3.3 percent at the design point inlet 
corrected flow rate (shown as a filled symbol) and over-predicted choking flow rate by 1.9 percent. The 
computation was run with ambient inlet pressure of 101,325 Pa, whereas the data was taken at 
approximately 86,000 Pa. The computation was performed with the measured clearances, which are 
minimum tip clearances. If clearances are highly non-axisymmetric in the test article, modeling the 
minimum tip clearance may explain part of the over-prediction in the computed performance. The 
questions of clearance sensitivity and Reynolds number sensitivity are left for future work. Instead, it was 
hypothesized that the over-prediction of performance may be due to inaccurate simulation of the 
unsteady, non-uniform impeller discharge flow and the subsequent impact on associated time-mean 
aerodynamic blockage at the diffuser throat. 

The scope of the research effort was narrowed to more closely examine the impeller exit flow. The 
focus of this paper is the CC3 impeller in its vaneless diffuser configuration. The compressor was 
modeled as a steady, periodic, single blade row computation in TURBO, which reduced computational 
expense and run-times. While the vaneless configuration obviates the impact of periodic unsteadiness, it 
was anticipated that lessons learned in modeling this simpler case may be useful in future work involving 
impellers and vaned diffuser configurations. The details of the compressor experiments are provided in 
the next section, followed by details of the computational effort. 
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Details of the Experiment 

The focus of the current work was the CC3 impeller in its vaneless diffuser configuration. A 
cross-section of the compressor rig is shown in Figure 3. Table 1 summarizes important design 
parameters. The impeller consists of 15 main blades and 15 splitter blades having 50° of backsweep from 
the radial direction. The design corrected speed was 21,789 rpm, and the measured minimum clearances 
at this speed were 0. 1524 mm (0.006 in.) at the leading edge, 0.6096 mm (0.024 in.) at the midchord 
“knee” and 0.2032 mm (0.008 in.) at the trailing edge (Ref. 7). At this clearance schedule and in the cold 
shape, the impeller hub at the trailing edge sat axially forward of the diffuser hub by 0.127±0.0508 mm 
(0.005±0.002 in.). This step in the hub was neglected in the computations. 

Total pressure and total temperature exit rakes were located at eight evenly spaced circumferential 
locations at radius ratio R/R TE = 1.18. Each total pressure rake had four elements located at 19.2, 43.4, 
68.2, and 93.1 percent of span, and each total temperature rake had three elements located at 23.1, 55.9, 
and 88.9 percent of span. Performance parameters were calculated using arithmetic averages of these rake 
elements. An inlet probe was used to acquire spanwise total pressure at a station 2.54 cm (1.0 in.) 
upstream of the impeller main blade leading edge (Fig. 4). Spanwise total pressure surveys were also 
acquired by using a constant-blockage probe immersed downstream of the impeller trailing edge at radius 
ratio R/R te = 1.10. Features which were neglected in the computations include a row of three inlet struts 
upstream of the impeller, impeller blade fillets, and backplate leakage paths. A set of 16 equally spaced 
5/ 16 th in. diameter pins, which held together the diffuser hub and shroud, were located in the test 
compressor at radius ratio R/R TE = 1.45. These structural elements were thought to be sufficiently far 
downstream of the exit rating station and were neglected in the computational model. 



Figure 3. — Cross-sectional drawing of the CC3 rig in the 
vaneless diffuser configuration (Ref. 7). A is the exit 
rating station, R/R TE = 1.18; 6 is the impeller main 
blade leading edge; C is the impeller splitter blade 
leading edge; D is the leading edge of a row of support 
struts; £ is the impeller trailing edge. 



(1.0 in.) upstream of impeller leading edge. 


TABLE 1.— CC3 IMPELLER DESIGN PARAMETERS 


Impeller main blade count 

15 

Impeller splitter blade count 

15 

Design corrected speed 

21,789 rpm 

Design inlet corrected flow rate 

4.54 kg/s 

Design exit corrected flow rate 

1.36 kg/s 

Impeller exit tip speed 

492 m/s 

Inlet blade height 

64 mm 

Exit blade height 

17 mm 

Leading edge tip radius 

105 mm 

Trailing edge radius, R TE 

215.5 mm 

Impeller backsweep 

50° 
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Details of the Computations 

The default turbulence model in TURBO Parallel Version 4 is the NASA/CMOTT k-s model with 
wall functions (Ref. 13). In addition to the default k-s model, Wilcox’s k-co turbulence model was used 
(Ref. 16). The implementation of the k-co model does not include wall functions, and the k-co model 
required a grid with sufficiently low y + so as to allow integration to the wall. Thus, two different grids 
were used in the current effort. A “baseline” grid with approximately 6.8 million points and y + on the 
order of 50 was used for k-e computations with wall functions enabled. A “refined” grid containing 
approximately 14 million points was also generated, allowing integration to the wall for k-co as well as k-e 
computations, in an effort to characterize the effect of varying grid resolution on aerodynamic 
performance. A histogram of y + for the k-s and k-co computations on the refined grid is shown in Figure 5. 
Typically, y + was less than 3 for a majority of the walls. Values of y + approached 5 along the blunt 
trailing edge of the impeller and at the shroud near the trailing edge of the suction side of the impeller, but 
these were within the viscous sublayer region. In all computational cases, inlet turbulence intensity was 
set to 1 percent, and inlet turbulent eddy viscosity ratio was set to 100.0. These values correspond to an 
inlet turbulence length scale of approximately 0.6 percent of the impeller inlet blade height. 

TURBO solves the unsteady RANS equations on multiblock structured grids. All grids, boundary 
condition files, and block connectivity files for TURBO were generated using the Turbomachinery 
Gridding System (TGS) (Ref. 17). The measured tip clearances, which were minimum tip clearances, 
were included in the computational grid. Impeller blades and tip clearances were gridded with O-blocks, 
and H-blocks were used to connect the O-blocks around the impeller main and splitter blades. The block 
topology of the refined grid is shown in Figure 6. The computational domain consists of the axial inlet, an 
impeller and splitter passage, the vaneless diffuser, and the S-bend to the axial outlet. A meridional view 
of the computed flow path is shown in Figure 7. 

Steady, single passage periodic simulations were converged at design speed using local time stepping. 
Solid walls were modeled with no-slip boundary conditions and were assumed to be adiabatic. The 
specific heat coefficient ratio was assumed constant. The inlet boundary condition was specified with 
spanwise profiles of total pressure, total temperature, radial flow angle, and tangential flow angle. The 
total pressure profile, shown in Figure 4, was the measured inlet probe survey data from the CC3 
experiment. The computations were referenced to an inlet pressure of 101,325 Pa, whereas the data was 
taken at reduced inlet pressures on the order of 85,500 Pa. Inlet total temperature was assumed constant 
288. 15 K, and zero inlet tangential flow angle was assumed. The profile of inlet radial flow angle was 
found by linear interpolation between the geometry angles of the hub and shroud at the inlet. The 
computations were throttled from choke, up the speed line, but were not taken to the stall boundary. The 
operating point was set by using a radial equilibrium exit boundary condition, specifying exit static 
pressure at the shroud. 



computations on the refined grid near design point. 
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Figure 6. — Block topology of the CC3 impeller and 
vaneless diffuser. 



Figure 7. — Meridional view of the CC3 impeller and 
vaneless diffuser model. 


Results 

Design speed computations of the impeller and vaneless diffuser were generated using TURBO. 

Three computed speed lines are compared to data in Figure 8. The three computed results are the k-s 
result with wall functions on the baseline grid, and the k-s and k-co results integrating to the wall on the 
refined grid. The computations were not taken to the stall boundary or deep into choke, as comparisons to 
the near-peak-efficiency data point were of interest in the current effort. The computed results were area- 
averaged in the circumferential direction at the exit rating station location, R/R TE = 1.18. The area- 
averaging of the steady computations in the rotating frame of reference is analogous to the time- averaging 
done by the steady rake instrumentation, which is fixed in the stationary frame of reference while 
immersed in the unsteady rotor exit flow. The circumferential averaging yielded a spanwise profile of 
computed total pressure and total temperature. An arithmetic average of the values of total pressure and 
total temperature at those spanwise locations which correspond to the rake element locations was 
calculated, which was analogous to the experimental averaging approach. 

The choking flow rates of the baseline and refined k-s computations agree to within 0.33 percent of 
each other. The k-co computation choked at corrected flow 0.5 percent lower than the refined k-s 
computation. This difference is attributed to a slightly larger growth in the boundary layer in the k-co 
result along the pressure side of the main impeller blade. The larger boundary layer results in increased 
aerodynamic blockage, reducing both effective flow area and choking flow rate. The k-co computation 
under-predicts choking flow rate by 1 .22 percent as compared to data. 

Both k-s results under-predict total pressure ratio and efficiency as compared to the data. The k-co 
computations better match the data on these performance parameters. The near-peak-efficiency data point 
and the computed operating points at nearby inlet corrected flow rates are shown as filled symbols in 
Figure 8. These operating points were compared. All three sets of computations over-predicted the work 
done by the impeller, as shown by the total temperature curves which are approximately 0.6 percent 
higher than the data. Compared to the near-peak-efficiency data point, the baseline k-s result was 

7.2 points lower in efficiency and 7 percent lower in total pressure ratio. The refined k-s result was 

8.3 points lower in efficiency and 8 percent lower in total pressure ratio than the data. The k-co result was 
1.2 points lower in efficiency and 0.6 percent higher in total pressure ratio than the data. In this case, the 
choice of turbulence model evidently had a much greater impact on the flow field than the selected level 
of grid density. 
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Figure 8. — Total pressure ratio, total temperature ratio, and adiabatic efficiency of CC3 vaneless diffuser configuration 
at design speed. Note: the computed speed lines were not run to stall. 


A comparison of shroud static pressures at these operating points is shown in Figure 9. The data have 
been normalized by the plenum pressure of 85,360 Pa and the computations are referenced to inlet total 
pressure 101,325 Pa. The shroud static pressures of the data, the k-co case, and both k-e cases show very 
good agreement from the inlet to the impeller trailing edge. Downstream of the trailing edge, the baseline 
and refined k-e computations have virtually identical shroud static pressures. The k-e computations have 
reduced pressure recovery through the diffuser as compared to the data and to the k-co computation. In the 
vaneless diffuser, the k-co simulation matches the data well through the initial half of the diffuser, but 
over-predicts pressure recovery towards the S-bend and beyond. 
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Figure 9. — Computed and experimental static pressure ratio along the shroud surface. Vertical lines mark geometry 
features of interest. 


In Figure 10 spanwise profiles of the computed flow angles (from tangential) and total pressures are 
compared against the profiles measured by the constant-blockage probe at radius ratio R/R TE = 1.10. This 
is a location nearly midway between the impeller trailing edge and the exit rating station, where the 
shroud static pressures of the k-co case and the k-e cases are diverging. The uncertainty of the measured 
flow angle is not obvious, so measured flow angles are shown with both ±1° and ±3° error bars for 
reference. The comparison shows that the baseline and refined k-s computations both yield nearly 
identical total pressure and tangential flow angle profiles. The flow angle profiles indicate that the k-e 
computations have a region of reversed flow (i.e., negative tangential flow angle) in the top 25 percent of 
the passage. This large region of reversed flow is an aerodynamic blockage which causes higher velocities 
and thus higher flow angles in the bottom 75 percent of the passage, as compared to the data. The 
reversed flow near the shroud is evidently the cause of the large deficit in total pressure observed in the 
k-s computations, which disagree with the measured pressure profile above 55 percent span. The data 
indicates reversed flow within the top 1 0 to 15 percent of span near the shroud. The flow angle profile 
from the k-co computation agrees with the data in this regard. The k-co computation over-predicts flow 
angle in the bottom 50 percent of span, but to less of a degree than either k-s computation. The k-co 
computation misses the pressure deficit seen in the survey data above 75 percent span. 

In Figure 1 1, a comparison of flow angles is made between the three computations at two locations; 
just downstream of the trailing edge (R/R T e = 1.01), where the shroud static pressures are in agreement, 
and at the exit rating station (R/R T e = 1-18), where there is a large difference in shroud static pressures. 
The flow angle profiles of the three computations match well near the trailing edge. At R/R T e = 1-01, the 
flow angle profiles of the k-e cases match each other to within ±0.5°, and the k-co case matches the k-e 
cases within ±1.5°, with slight disagreement above 95 percent span. The k-e cases indicate reversed flow 
beginning at 97 percent span, the k-co case indicates reversed flow at 95 percent span. There is large 
disparity between the k-e cases and the k-co case at radius ratio R/R te = 1.18. The flow angle profiles of 
the two k-e computations, which agree to within ±1.0°, indicate that the flow has reversed in nearly 
30 percent of the span near the shroud. The k-co case has reversed flow in only the upper 5 percent of span 
near the shroud. 
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Figure 10. — Area-averaged computed total pressure ratio (left) and flow angle (right) compared to survey data at 
radius ratio R/Rje = 1.10. 



Figure 1 1 . — Area-averaged flow angles comparing results of the k-£ and k-u) turbulence models at radius ratios 
R/Rte = 1.01 (left) and R/R TE =1.18 (right). 
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Figure 12. — Radial velocity contours (normalized by 287.6 m/s) at the trailing edge of the impeller, R/Rte = 1 .01 (left), 
and at the exit rating station R/Rte =1.18 (right). Impeller blade rotation is right-to-left. 

The deficit in static pressure recovery seen in the k-s computations is attributed to the large region of 
reversed flow. This region is shown in Figure 12, which plots the computed radial velocity contours near 
the impeller trailing edge, R/R T e = 1.01, and at the exit rating station, R/R TE = 1.18. The contours near the 
trailing edge show the splitter blade in the middle of the plot, and the main blade at the left and right sides 
of the plot. At R/R te = 1.01, the low-momentum regions near the hub that extend into mid-span are 
associated with the impeller blade wakes. The low radial velocity region at the shroud is the tip clearance 
flow. The spanwise extent of the clearance flow is slightly larger in the k-s computations as compared to 
the k-co computation, indicating a stronger tip clearance vortex in the k-s computations. The flow fields of 
the two different turbulence models are qualitatively and quantitatively similar at the trailing edge. If 
performance is computed from the inlet plane to the plane defined by R/R TE = 1.01, the refined k-s and 
k-co computations agree to within 0.7 percent in total pressure ratio and to within 1.0 point in adiabatic 
efficiency. At the exit rating station R/R TE = 1.18, the two k-s computations are qualitatively very 
different from the k-co computation, and the computed performance parameters at R/R TE =1.18 amount to 
a 8.4 percent difference in total pressure ratio and 7.0 point difference in adiabatic efficiency between the 
refined k-s and k-co computations. 

Figure 13 plots entropy on planes of constant radius through the last 30 percent of impeller chord to 
show loss development through the passage. In the figure, the pressure side of the main blade is on the 
right, the splitter blade is in the middle, and the suction side of the main blade is on the left. The 
computations differ at the comer formed by the pressure side of the main blade and the shroud. The 
higher magnitude of entropy in the refined k-s case indicates larger loss due to a stronger tip clearance 
flow than in the baseline k-s case. The spanwise extent of the high entropy region is smaller in the k-co 
case, indicating a weaker tip clearance flow than in the k-s cases. Qualitatively, the entropy contours 
among the three computations are similar through the trailing edge. 
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Figure 13. — Entropy contours at constant-radius planes over the last 30 percent of impeller chord, with the splitter 
blade shown in the center of the passage. 
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Figure 14. — Area-averaged velocity contours (normalized by 
287.6 m/s) showing streamwise components (radial component 
in the radial diffuser and axial component after the radial-to- 
axial S-bend). 

Figure 14 plots radial velocity through the radial diffuser and axial velocity after the S-bend. These 
velocities were area-averaged in the circumferential direction for the k-s and k-co results on the refined 
grid. This shows the extent of the reversed flow through the diffuser along the shroud in each case. In 
both computations, the impeller exit flow is similar, but differences arise in the diffuser. Both computa- 
tions show reversed flow at the shroud at close proximity to the impeller trailing edge. In the k-s 
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computation, the reversed flow persists up the radius of the diffuser along the shroud surface. In the k-co 
case, there is a relatively small bubble of reversed flow which weakens as radius increases through the 
diffuser. The separated region in the k-e cases reattaches through the S-bend and there is no negative axial 
velocity at the end of the computational domain. The diffuser flow field had large sensitivity upon the 
choice of turbulence model. A discussion on the differences in the diffuser flow field is presented next, as 
well as comments and suggestions for future work. 

Discussion 

A better understanding of the cause of the disparate diffuser flow fields between the k-e and k-co 
computations is of interest for future work. The flow fields indicate that the k-co computation has greater 
turbulent mixing in the diffuser, which causes the shroud separation to weaken and mix out in a manner 
that is not occurring in the k-e computations. Figure 15 shows a band of high turbulence kinetic energy 
(TKE) in the diffuser of the k-e case which is not evident in the k-co case. This region is attributed to the 
shear layer at approximately 70 percent span due to the interface of the separated flow and core flow in 
the k-e case. Generally, it was observed in the k-co computation that the levels of both production and 
dissipation of TKE were higher than the k-e computations near the shroud within the impeller passage, 
but these terms were lower than the k-e computations through the diffuser. 

It may be of interest in future work to generate a simulation using a standard k-e model, which sets 
constant C M = 0.09 in the turbulent eddy viscosity equation rather than a computed value as in the 
CMOTT k-e model. As shown in Figure 16, the computed C M in the refined CMOTT k-e result yields 
values lower than 0.09 in the shear layer region at 70 percent span in the diffuser at the exit rating station. 
A standard k-e model would increase C M in this region, resulting in increased turbulent eddy viscosity and 
turbulent mixing, which may mix out or weaken the separation observed in the current k-e cases. 



Figure 15. — Turbulence kinetic energy at constant-radius planes from R/Rte : 
is shown in the center of the passage. 


0.9 to R/Rte = 1.18. The splitter blade 
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Refined k-epsilon, R/R_TE=1.18 



Figure 16. — C M of the refined grid k-£ computation at 
radius ratio R/Rte =1.18 (exit rating station). 


The current results suggest that the flow developing through the radial diffuser may largely depend on 
the selection of turbulence model. In a vaned diffuser configuration, a strong separation entering the 
diffuser impacts the aerodynamic blockage, which affects the pressure rise and flow capacity in the 
diffuser. The impact of different turbulence models on this separation, and thus on the stage performance, 
should be analyzed in future work. The vaned diffuser computation poses a more complicated problem, as 
the blockage will also be dependent upon the unsteadiness due to the impeller/vane interaction in a vaned 
diffuser configuration. 


Conclusion 

A computational study of the CC3 impeller in its vaneless diffuser configuration was undertaken. 
TURBO CFD simulations were generated on: (1) a baseline grid using the default CMOTT k-s turbulence 
model, (2) a refined grid using the default CMOTT k-s turbulence model, and (3) a refined grid using 
Wilcox’s k-co turbulence model. The results indicated that all three computational cases matched the data 
well from the inlet, through the impeller, and up to the trailing edge. The area-averaged flow angle 
profiles of the baseline and refined k-s cases matched to within ±0.5° just downstream of the impeller 
trailing edge (R/Rte = 1.01). At this radius ratio, the flow angle profile of the k-co case matched both k-s 
cases to within ±1.5°, and the refined k-s and k-co computations agreed to within 0.7 percent in total 
pressure ratio and to within 1.0 point in adiabatic efficiency. Beyond the trailing edge and into the 
diffuser, the k-co computation continued to match the data, but the k-s computations yielded pressure and 
efficiency deficits as compared to data due to a large reversed flow region extending from the shroud 
surface of the vaneless diffuser. At the experimental exit rating station (R/Rte = 1 . 1 8), the baseline k-s 
computation is 7.2 points lower in efficiency and 7 percent lower in total pressure ratio than the data. The 
refined k-s computation is 8.3 points lower in efficiency and 8 percent lower in total pressure ratio than 
the data. The k-co computation is 1.2 points lower in efficiency and 0.6 percent higher in total pressure 
ratio than the data. It was shown that the potential impact of the choice of turbulence model may be quite 
large on the diffuser flow field. 
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